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We investigate the Landau parameters for the instabilities in spin and charge channels in the 
nondegenerate extended Hubbard model with intersite Coulomb and exchange interactions. To this 
aim we use the spin rotationally invariant slave boson approach and we determine the necessary 
inverse propagator matrix. The analytically derived spin Landau parameter Fq for the half filled 
band uncovers the intrinsic instability of the nondegenerate Hubbard model towards ferromagnetism 
— negative intersite exchange interaction triggers a ferromagnetic instability at half filling before the 
metal-insulator transition, indicated by the divergence of the magnetic susceptibility at Fq = — 1. 
This result is general and the instability occurs in the strongly correlated metallic regime for any 
lattice, in three or two dimensions. Next as an illustrative example we present numerical results 
obtained for the cubic lattice with nearest neighbor exchange J and Coulomb V elements and 
arbitrary electron density. One finds that the range of small doping near half filling is the most 
unstable one towards spin polarization, but only in the case of ferromagnetic intersite exchange 
J < 0. Charge Landau parameter Fq is lowered near half filling by increasing U when the intersite 
Coulomb interaction V is attractive, but in contrast to Fq at J < 0 it requires an attraction beyond 
a critical value V c to generate the divergence of the charge susceptibility at Fq = — 1 in the metallic 
phase. This instability was found for a broad range of electronic Filing away from half filling for 
moderate attraction. 


I. INTRODUCTION 

The Hubbard model has been initially introduced, in¬ 
ter alia , to describe metallic magnetism 00 - This is 
frequently forgotten in recent research which focuses on 
other instabilities which are hidden and not so easy to in¬ 
vestigate, as the instabilities of the Fermi surface or the 
onset of superconductivity @. The ferromagnetic (FM) 
state is easily derived from the Stoner criterion and oc¬ 
curs at sufficiently large local Coulomb interaction U. 
The Stoner instability of a nonmagnetic state is deter¬ 
mined using the Hartree-Fock approximation and is thus 
predicted for any lattice and filling of a nondegenerate 
band. However, it was soon realized that it does not 
imply ferromagnetism as electron correlations strongly 
renormalize the Stoner parameter which follows from the 
Hartree-Fock approach and remove the FM instability 
in most cases. The metallic state of strongly correlated 
electrons in the nondegenerate Hubbard model is a cor¬ 
related Fermi liquid which is characterized by the Fermi 
liquid parameters. In the frequently considered model 
of a rectangular density of states, the FM instability is 
absent. This conclusion follows also by considering the 
Fermi-liquid antisymmetric parameter Fq, which lower 
bound was found to be —3/4 0. Indeed, ferromagnetism 
in a nondegenerate band is rare and only very few sys¬ 
tems exist in nature which could be classified as itinerant 
ferromagnets with a single band, such as a weak itinerant 
ferromagnet ZrZn 2 Q. 

After several decades of research on the Hubbard 
model, it became clear that the mechanism of itinerant 


ferromagnetism is subtle and other effects are needed to 
stabilize it. One of them is high degeneracy of the ground 
state Q realized when the band is flat 0-0. Ferromag¬ 
netism in realistic flat-band systems has been proposed 
for a triangular lattice of Naa,Co 02 EH and for the p or¬ 
bitals in honeycomb lattices with ultracold atoms E3. In 
a two-dimensional (2D) Hubbard model on a square lat¬ 
tice FM instability occurs at the van Hove singularity of 
the density of states [l^-[T3|. It becomes relevant away 
from half filling when the band is flat due to the next 
neighbor hopping t'. This follows from a more general ap¬ 
proach which predicts the FM instability by investigating 
the divergence of the uniform magnetic susceptibility. 

Other factors which stabilize itinerant ferromagnetism 
in real materials go beyond the Hubbard model itself 
these are: (i) direct intersite exchange coupling jl5l - fl7j . 
or (ii) band degeneracy with active Hund’s exchange cou¬ 
pling EMI. The second mechanism is better known 
it is responsible for ferromagnetism in transition metals 
Fe, Co and Ni which have rather large magnetic mo¬ 
ments and FM states occur in a degenerate 3d band 
where Hund’s exchange plays a prominent role 0] . Also 
here the Stoner parameter is strongly renormalized by 
the electron correlation effects E0> but the FM insta¬ 
bility survives and is much easier to obtain than for the 
nondegenerate s band EH - FM states have been investi¬ 
gated for realistic transition metals using different meth¬ 
ods including more recent studies with the Gutzwiller 
wave function (2lj . This subject is interdisciplinary and 
itinerant ferromagnetism plays also a role in a ferro- 
magnet/insulator/superconductor ballistic junction |22| . 
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while the essential role played in the onset of ferromag¬ 
netism by the orbital degeneracy was also demonstrated 
recently for the multiband Hubbard models on square 
and cubic lattices [23] relevant both for p-orbital bands 
with ultracold fermions in optical lattices, and for elec¬ 
tronic 3d-orbital bands in transition-metal oxides [24]. 

Here we focus on the extended nondegenerate Hub¬ 
bard model which describes the simplest correlated met¬ 
als and investigate the influence of intersite Coulomb and 
exchange interactions on the charge and magnetic insta¬ 
bilities. Thereby we follow the route initiated by Voll- 
hardt [4] who investigated Landau Fermi liquid parame¬ 
ters Fq and Fq for 3 He. Further studies gave the Fermi 
liquid interaction and the quasiparticle scattering ampli¬ 
tude on the Fermi surface in 3 He [25[. The perturbative 
method predicts the change of sign of the parameter Fq 
in the 2D square lattice for small U ~ 2f [26[, where 
t is the hopping parameter. The calculation within the 
slave boson method in the doped 2D Hubbard model [27] 
does not confirm this result and gives a negative Fq in 
the entire range of U and its saturation with increas¬ 
ing U towards a negative value but larger than —1 and 
nearly independent of the doping. This behavior con¬ 
firms the earlier result that the nondegenerate Hubbard 
model does not exhibit a FM instability even when U is 
very large [4] . We remark that saturated ferromagnetism 
near half filling and at U —> 00 which stems from the 
Nagaoka state is nevertheless not excluded, as shown by 
several studies @, |28[ including a slave boson approach 
[29l ]. This subject goes however beyond the scope of the 
present study. 

Since Mott insulating ground states arise at large U 
and at half filling, we perform our investigations in a 
framework which is able to capture interaction effects 
beyond the physics of Slater determinants, ft is an ex¬ 
tension of the Kotliar and Ruckenstein slave boson repre¬ 
sentation that reproduced the Gutzwillcr approximation 
on the saddle-point level (30]. It entails the interaction 
driven Brinkman-Rice metal-to-insulator transition |3lj . 
A whole range of valuable results have been obtained with 
Kotliar and Ruckenstein [30| and related [H,[33[ slave bo¬ 
son representations which motivate the present study. In 
particular they have been used to describe antiferromag¬ 
neti c Irl . ferromagnetic [13, spiral [35l - l37j and striped 
[38l - l4l| phases. Furthermore, the competition between 
the latter two has been addressed as well [HI] . The in¬ 
fluence of the lattice geometry on the metal-to-insulator 
transition was discussed, too [43|. For instance, a very 
good agreement with Quantum Monte Carlo simulations 
on the location of the metal-to-insulator transition for the 
honeycomb lattice has been demonstrated [3]. Finally, 
further motivation comes from the strongly inhomoge¬ 
neous polaronic states that have been found in correlated 
heterostructures using the Hubbard model extended with 
intersite Coulomb interactions [45| . 

Furthermore, comparison of ground state energies to 
existing numerical solutions have been carried out for the 
square lattice, too. For instance, for U = 4f it could be 


shown that the slave boson ground state energy is larger 
than its counterpart by less than 3% [35| . For larger 
values of U, it has been obtained that the slave boson 
ground state energy exceeds the exact diagonalization 
data by less than 4% (7%) for U = 8t (20 1) and doping 
larger than 15%. The discrepancy increases when the 
doping is lowered ■ It should also be emphasized that 
quantitative agreement to quantum Monte Carlo charge 
structure factors was established [46| . 

The purpose of this paper is to derive and evaluate 
Fermi liquid Landau parameters Fq and Fq for the metal¬ 
lic state in the extended Hubbard model. Despite of the 
above discussion on the FM instabilities, the influence of 
intersite Coulomb and exchange interactions on Landau 
parameters was not analyzed in the Hubbard model un¬ 
til now. We investigate the instabilities towards uniform 
FM order and charge instabilities using the spin rotation 
invariant (SRI) representation [HIHH of the Kotliar and 
Ruckenstein slave boson approach. We show analytically 
that weak FM intersite interactions are sufficient to trig¬ 
ger a divergence in the magnetic susceptibility expressed 
by the antisymmetric Landau parameter Fq which im¬ 
plies that the FM instability occurs due to such interac¬ 
tions for any lattice. For a representative example of a 
three-dimensional (3D) cubic lattice we present also nu¬ 
merical analysis, and compare them at half filling with a 
two-dimensional (2D) square lattice. 

The paper is organized as follows: The extended Hub¬ 
bard model is introduced in Sec. El together with its 
Kotliar and Ruckenstein SRI slave boson representation. 
In Sec. uni we perform the saddle-point approximation 
and present the resulting system of coupled nonlinear 
equations. Fluctuations are captured within the one-loop 
approximation described in Sec. IIV1 which allows one to 
determine analytically Fermi liquid Landau parameter 
Fq at half filling. Numerical results are presented and 
discussed in Sec. El where we address first Fq and the 
instabilities towards uniform spin polarization (ferromag¬ 
netism) and next consider Fq and charge instabilities in 
Secs. IV Al and IVBl The paper is summarized in Sec. ED 


II. EXTENDED HUBBARD MODEL 

Numerous studies of correlated electrons have been de¬ 
voted to the properties of the Hubbard model on a square 
lattice, especially after Anderson’s proposal that it rep¬ 
resents a minimal model for the d electrons within the 
Cu0 2 layers common to the high T c superconductors 
[47l ] . Yet the Hubbard model assumes a perfect screen¬ 
ing of the long-range part of the Coulomb interaction. 
This may be questionable and the relevance of this ap¬ 
proximation may be assessed by considering the extended 
Hubbard model that reads: 

H ~ ^' i .i C ia C jo- + \ ^ ~ 2 / 

i,j,cr i ' ' ^ ' 
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2 5Z n 0(i n j )+ ^ (i) 


and includes intersite Coulomb 14, and exchange Jij in¬ 
teractions. These elements decay fast with increasing 
distance | Ri — Rj |, but extend in general beyond nearest 
neighbors. Here c\ a are electron creation operators at site 
i with spin a, m a = c' icT c i(T . and 5, are spin operators. 
We consistently use the particle-hole symmetric form for 
both density-density interaction terms. 

Although one expects that V % j >0, in certain cases ef¬ 
fective intersite Coulomb interactions may be attractive 
[481 ]. Therefore, we shall treat {Vjj } as effective parame¬ 
ters and consider both signs of them below. Similar, for 
the exchange elements {Jij} we shall also consider both 
antiferromagnetic (Jij > 0) and ferromagnetic (Jij < 0) 
exchange elements. In a transition-metal oxide the dom¬ 
inating interaction is the superexchange and the antifer¬ 
romagnetic coupling is expected. On the contrary, in a 
metallic system the direct exchange elements which arise 
from the Coulomb interactions could be ferromagnetic. 
Thus, the present theory makes predictions concerning 
itinerant magnetism. 

In the original Kotliar and Ruckenstein representa¬ 
tion the spin interaction term remains a four-fermion 
term. This is no longer the case in the SRI represen¬ 
tation [32|, HIJ which we therefore adopt and extend for 
our study. In this framework the Hamiltonian Eq. © 
may be represented as: 


H = Y Y 

i,j,cr gg' g\ 

+ u'£(d\d l -Wflf!<,' + ] 
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1 + * i-E&zJ 


E4- y ; ^ [ T PP'PjppiPjpip' ’ 


i,j gg'gi 


( 2 ) 


where we introduced: 


z = e^L M R p + p^L M R d, (3) 


with P = \ 0 PpJ^‘ an< 4 1f being the Pauli matrices. 

We further use p a(J , = cra'p_ cr , _ cr , and 


M = 


1 

2 


1 +e ] e + '^p\p ll +d ] d 




L = [(l — d^d) 1 — 2p^p ] 2 , 
R= [(l — e f e) 1 — 2ptp ]~^ . 


(4) 


For more details see Ref. [49(. In Eq. © we used the 
mapping of the fermionic degrees of freedom onto bosons 


and expressed the hole doping operator as follows, 


Y t = e| e i - d\ d t , 

(5) 

and the spin operator as 


^ ^ v ^^'PlGGlPiGlG' ' 

(6) 


gg'gi 


Since the auxiliary operators span an augmented Fock 
space, physically meaningful results may be obtained pro¬ 
vided they satisfy local constraints. For each site they 
read: 

e t e + ^p^ + d t d= 1, (7) 

+ 2 dU, (8) 

G [I 

Y ft 1 7°fa = POP + P*P0 ~ * P ■ (9) 

G,G r 

They may be enforced in path integral formalism. 

We remark that any slave boson representation pos¬ 
sesses an internal gauge symmetry group. In our case it 
may be made use of to simplify the problem and to gauge 
away the phases of the e and p M bosons by promoting 
all constraint parameters to fields [33| . leaving us with 
radial slave bosons fields [50]. Their approximate values 
that are obtained in the saddle-point approximation may 
be viewed as an approximation to their exact expectation 
values that are generically non-vanishing a The slave 
boson field corresponding to double occupancy d has to 
remain complex however as emphasized by several au¬ 
thors maa 

Besides, the saddle-point approximation is exact in the 
large degeneracy limit, and the Gaussian fluctuations 
provide the 1/N corrections [Hj. Moreover it obeys a 
variational principle in the limit of large spatial dimen¬ 
sions where the Gutzwiller approximation becomes exact 
for the Gutzwiller wave function [54j . Furthermore, it 
could be shown in this limit that longer ranged inter¬ 
actions are not dynamical and reduce to their Hartree 
approximation [551 ]. Therefore, our approach also obeys 
a variational principle in this limit when applied to the 
above extended Hubbard model Eq. ©. 

One can infer from these formal properties that the 
approach captures characteristic features of strongly cor¬ 
related electrons as the suppression of the quasiparticle 
residue and the Mott-Hubbard/Brinkman-Rice transi¬ 
tion [3lJ to an insulating state at half filling with increas¬ 
ing on-site Coulomb interaction. Of particular relevance 
is the influence of the longer-ranged Coulomb interaction 
on the latter transition. 


III. SADDLE-POINT APPROXIMATION 

Though such functional integrals can be calculated ex¬ 
actly for the Ising chain [50:] and some toy models [5lj . 






even with the Kotliar and Ruckenstein representation 
[5(1 ]. this is unpractical on higher dimensional lattice. 
Here we rather resort to the saddle-point approximation. 
In the translational invariant paramagnetic phase all the 
local quantities are site independent, and the action at 
saddle-point reads (/3 = 1/ksT), 

S = (3L (s B + S F + , (10) 

with 

S B = a(e 2 + d 2 + p 2 0 - 1) - /3 0 (Po + 2rf2 ) 

+ Ud 2 + ^V 0 Y, (11) 

S F = -- p Y j \n(l + e~^) . (12) 

k,a 

For the extended Hubbard model 0 the quasiparticle 
dispersion in Eq. m reads: 

E ka = + A) -\u- \v 0 Y - p , (13) 

in which we introduced the Fourier transform of the in¬ 
tersite Coulomb repulsion, 

v k = jJ2 v v e ~ it{ * j ~ Rz) - ( 14 ) 


It is convenient to define the Fourier transform of the 
intersite exchange elements in a similar way, 

= (!5) 


From Eq. m one obtains the following set of saddle- 
point equations: 

Pq + e 2 + d 2 — 1 = 0, 
pi + 2d 2 = n, 

1 dz% 1 , .1 dY , 

^5+jW(16) 

^-p-e = p 0 -a, 

2p 0 dpo 

i f ) z 2 1 1 r)Y 

= ! ^-a) + c- u . 

Here we have introduced the averaged kinetic energy, 


£ = 


J dujp{uj)ujf F (zqU + fa - ^(7 - ivb Y - p 


(17) 

which determination involves the density of states p{co) 
and f F is the Fermi function. It turns out that the last 
three equations in the set Eqs. (fl6l) may be merged into 
a single one, 


y ( 1 - y) 

n( 1 — §)ed 


e + U = 0. 


(18) 



n 


FIG. 1. (Color online) Average kinetic energy per site e for 
the extended Hubbard model as obtained for the 2D square 
lattice (blue, solid line) and for the 3D cubic lattice (red, 
dashed line) as functions of the electron filling n. Parameters: 
t= 1. 


Here we have introduced 

y=(e + d) 2 , (19) 

in terms of which the saddle-point equation finally reads: 

V 3 + (u - 1 )y 2 = uS 2 , (20) 

where S = 1 — n is the doping away from half filling 
and u = U/Uo is the Coulomb parameter in units of Uq 
defined through 


= ( 21 ) 

At half filling a metal-to-insulator transition occurs for a 
given lattice at U c , which is defined as follows 

U c = lim U 0 = — 8e. (22) 

S-s-0 

It is very remarkable that neither Jy nor Vij elements 
enter Eq. (1201) . Hence, in a paramagnetic phase, the in¬ 
tersite interactions only influence the fluctuations and do 
not change electron localization due to strong onsite in¬ 
teraction U. In particular, the nearest neighbor Coulomb 
interaction V has no influence on the Mott gap as dis¬ 
cussed by Lavagna [57]. Furthermore the double occu¬ 
pancy is in exact agreement with the Gutzwiller approx¬ 
imation as derived by Vollhardt, Wolfle and Anderson 
fhH . In the present case of a 3D cubic lattice the double 
occupancy vanishes at half filling for U c = 16.0387t (at 
the Brinkman-Rice transition |3l|). 
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This behavior is generic and the transition occurs for 
other lattices in a qualitatively equivalent way. As an 
example we take also the case of a square lattice, where 
the number of nearest neighbors is reduced from z 1 = 6 
to z\ = 4. One expects that the value of U c (1221) will 
roughly scale with the number of nearest neighbors z i, 
and corrections follow from the shape of the density of 
states. In fact, the averaged kinetic energy |ej is reduced 
somewhat less than by a factor 2/3 from that found in 
the cubic lattice due to the shape of the 2D density of 
states which is closer to a rectangle than in the 3D case, 
see Fig. [T] In fact, for the 2D square lattice the metal- 
to-insulator transition occurs at U c = 2(8/7r ) 2 t. 

At the critical value of U = U c the effective mass m* 
diverges in the paramagnetic state (we recall that in the 
strongly correlated regime local moments are formed) 
and its inverse, 


being quasiparticle residue, vanishes, see Fig. [2] Here m 
is the electron mass and z 2 stands for the reduced jump 
in the electronic filling, (n% ), when the momentum k 
changes from inside to outside of the Fermi surface [ij. 
The plot shows that only for filling n > 0.8 the Fermi 
liquid is substantially renormalized in the regime of large 
U > U c . 

Solving the saddle-point equations (1161) yields at half 
filling, 



* 

m/m 



n 


FIG. 2. (Color online) Inverse effective mass Zg (1231) for the 
extended Hubbard model on the cubic lattice. 

Here we use the notation k = (k,w). Performing the 
calculation to one-loop order, one can make use of the 
propagator Sij(k ) as given in the Appendix. The suscep¬ 
tibilities are given by the matrix elements of the inverse 
propagator as follows: 

Xs(k) = 2p 0 S' 77 (fc), 

Xc(k) = 2 e 2 S^(k) - AedS^(k) + 2 d 2 S^(k). (28) 


The present analysis shows that its doping dependence is 
universal in the extended Hubbard model, and is influ¬ 
enced neither by nor by Jij. 


IV. ONE-LOOP APPROXIMATION TO THE 
RESPONSE FUNCTIONS 

Having mapped all degrees of freedom onto bosons al¬ 
lows us for directly evaluating the spin and charge re¬ 
sponse functions. Indeed, following Ref. 0 
and density fluctuations may be expressed as 


SN = Sric = S^d — e^e), 

a 

in the SRI representation. The spin and charge autocor¬ 
relation functions can be written in terms of the slave 
boson correlation functions as: 

Xs{k)= {5n a {-k)5n a > (k)} = (6S(-k)5S(k)), 

a,a 

Xc(k)=^(6n c ,(-k)6n ( ,'( k )) = (SN(-k)SN(k)). (27) 


the spin 

Xs{k) 

(25) 

where 

(26) 



As emphasized and analyzed by several authors, see e.g. 
Refs. [59L16011 • the Fermi liquid behavior is obtained when 
considering the above Xs{k) and Xc(k) in the long wave¬ 
length and low frequency limit. However, in contrast to 
the conventional random phase approximation (RPA) re¬ 
sults, the obtained Landau parameters involve effective 
interactions, which differ in the spin channel and in the 
charge channel. 

Explicitly the spin susceptibility reads: 

_ Xo(k) _ 

1 +A^ X o{k) +Bxi{k) +C[xl(k) ~Xo(k)x 2 (k)]' 

(29) 




B = 

C = 


a - /?o +££ (|^) l 




z 0 ( dz t 


Po \dp 3 

z 0 \ 2 f dz t 


Using 

2 

OP3 


2p 0 J \dp 3/ 

(e — d)( 1 — 5 2 ) — 2p^5(e + d) 
(1 - < 5 2 ) 3/2 


+ 7 J E , (30) 

( 31 ) 

(32) 


(33) 
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= (1^575 + ■*) [1 - ^ + 2J=)] 

- 4<5(e-d)(l-d 2 )} . (34) 

together with the saddle-point Eqs. (fTUl) yields the Lan¬ 
dau parameter Fq at half filling as: 

cm _ o^Ar(o) u ( l ~ u )( 2 + u) + (1 + u)(J 0 /8e) /otn 

- Z£I ^f ttttjo ?i it; > ia 0 / 


where we have introduced the bare density of states at 
the Fermi energy Np'* = p(Ep). We decompose Eq. (l35l> 
into a regular and singular part, and using Eq. (1221) we 
find, 


Fq = 2Np' > e 


u( 2 + u) 
(1 + u) 2 


Jo/Uc \ 

1 — u 2 J 


(36) 


The latter term depends solely on the k = 0 component 
of the exchange coupling, Jq = J^ =0 defined in Eq. m- 
We emphasize that the ferromagnetic instability deduced 
from Eq. (1361) is general and occurs in all cases below the 
metal-insulator transition when J 0 < 0. This remarkable 
result follows from the band narrowing when U —> U c 
which amplifies the effects of the intersite exchange in¬ 
teraction. 

Moreover, spin and charge fluctuations separate at the 
one-loop order, and the intersite Coulomb elements in 
(HU) have no effect on the value of Fq. Unfortunately, a 
similar analytic result for Fq could not be derived at half 
filling as the dependence on Vq = Vg_ 0 defined in Eq. 
m is not sufficiently transparent. 

Equation (l36l) is the central analytic result obtained 
at half filling. Its physical origin lies in the fact that, in 
the limit of vanishing hopping, the Hubbard model favors 
the formation of localized magnetic moments that order 
according to the exchange couplings, for instance ferro- 
magnetically for Jo < 0. Our result indicates that a min¬ 
imum of coherence of the quasi-particles z F is necessary 
to destabilize the ferromagnetic order. It only depends 
on j o = Jq/U and, for a rectangular density of states 
(DOS), reads: 


Jl 4jo + jl + (1 - jo)yJl - 6jo + Jo - 1 /oh7 ^ 

=- M -• (37) 


which behaves as z 2 F — — 2jo for small ferromagnetic ex¬ 
change. Hence, for Jo —> 0“ the ferromagnetic instability 
takes place at U = U c , while it is absent for J 0 = 0. As a 
result we have shown that introducing an arbitrarily weak 
exchange Jo ^ 0 interaction results in a singular behav¬ 
ior of Fq in the vicinity of the Brinkman-Rice point U c 
(u = 1). This is in marked contrast to the known results 
for a flat band in the absence of intersite exchange: 


F 0 a = -1 


1 


(1 + U/U c ) 2 


(38) 


While the simple form 


U(2U C -U) 1 

0 (U c - U) 2 + (1 - U/U c ) 2 


( 39 ) 


could be derived in absence of Coulomb elements this 
is no longer the case when they are taken into account. 
Note that the property Fq(U) = Fq(—U) can be derived 
on a more general ground [4j. 


V. NUMERICAL RESULTS 
A. Ferromagnetic instability — Fq parameter 

In this Section we compute the Landau parameter Fq 
for the extended Hubbard model with nearest neighbor 
exchange or Coulomb interactions on the cubic lattice 
to analyze its filling dependence. We first investigate 
half filling (n = 1), evaluating the obtained analytic ex¬ 
pression Eq. (l36l) to demonstrate the divergent behavior 
at finite intersite exchange. The Hubbard model does 
not exhibit the FM instability in the metallic regime as 
Fq > —1, see Fig. [3] However, as anticipated above, an 
infinitesimal FM coupling Jq generates an instability of 
the nonmagnetic state at half filling when the Coulomb 
interaction U approaches U c . 

Although Eq. (15C)1) is more general, we shall consider 
below the extended Hubbard model with nearest neigh¬ 
bor exchange interactions only, 

_ J J for a bond (■ ij ) with j e A f(i) , 

Jij ~\0 otherwise ’ ^ 

where Af(i) stands for the set of nearest neighbors of 
site i. We have found that the model is unstable at half 
filling against the FM order near the metal-to-insulator 
transition. The location of the instability depends rather 
sensitively on the FM coupling, from U c for J = 0 + 
down to 0.33[/ c f° r J/U = —0.2. Larger values of —J/U 
are nonphysical and they are not considered below. On 



FIG. 3. (Color online) Landau parameter Fq for the extended 
Hubbard model at half filling on the cubic lattice. Different 
lines from top to bottom for decreasing J/U = 0.2, 0.1, 0.04, 
0, -0.04, -0.1, -0.2. 
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FIG. 4. (Color online) Landau parameter Fq for the Hubbard 
model on the cubic lattice. The FM instability occurs only 
near half filling (at the boundary of the white region where 
Fq < — 1) at very large U > 32 1 (U > 2U C ). 


the contrary, a positive (AF) intersite exchange coupling 
Jo > 0 suppresses the tendency towards ferromagnetism, 
and the parameter Fq becomes positive and large near 
the metal-to-insulator transition, see Fig. [31 

Since the Landau parameter Fq is directly propor¬ 
tional to the product of the bare kinetic energy and the 



FIG. 5. (Color online) Instability lines of the unpolarized 
state towards FM order as given by the divergence of the 
magnetic susceptibility (Landau parameter Fq = —1) for the 
extended Hubbard model with FM exchange J < 0 on the 
cubic lattice. Different lines from top to bottom for: J/U = 0, 
-0.01, -0.05, -0.1, -0.15, -0.2. 
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FIG. 6. (Color online) Landau parameter Fq for the extended 
Hubbard model with AF J/U = 0.1 on the cubic lattice. 
At sufficiently large U the parameter Fq is positive (white 
region). 

density of states at Ep, eNp\ it shows a lesser sensitiv¬ 
ity to the band structure. For instance, one finds at half 
filling eNp' 1 = —0.5 for a flat (rectangular) DOS, and 
eNp' 1 ~ —0.57 for the cubic lattice. This gives the value 
of Fq ~ —0.86 at J = 0. We have verified that sim¬ 
ilar values are obtained for other typical DOSs at half 
filling, except for the 2D square lattice, where the van 
Hove singularity at n = 1 results in the FM instability at 
infinitesimal U = 0 + El - disregarding the stronger AF 
instability following from nesting. 

Next we consider the doping dependence of the Landau 
parameter Fq in the Hubbard model, i. e., at Jij = 0. 
As one can see in Fig. [4] Fq reaches the critical value 
Fq = — 1 for U — 32 1 at S = 0 + , namely in the regime 
where Nagaoka ferromagnetism is expected. The critical 
U increases linearly with doping. Furthermore, the cusp 
in the density of states for the cubic lattice at n ~ 0.43 
makes itself clearly visible in Fig. [H with reduced values 
of \Fq | at lower electron density. 

Also away from half filling, finite FM exchange cou¬ 
pling Jo < 0 triggers the FM instability Fq = — 1 at 
significantly lower values of U. For instance, already 
J/U = —0.01 brings this instability down to the values 
of U ~ 20 1 for the doping S < 0.57 (filling n > 0.43) 
where the DOS is almost independent of energy, see 
Fig. [5] When J/U = —0.05, the FM instability occurs at 
U < lOt in the same doping regime, and comes down also 
for the lower electron fillings. For lower J 0 the FM insta¬ 
bility occurs at even lower values of U. This is in contrast 
to the calculations to the two-band model, where the FM 






















FIG. 7. (Color online) Landau parameter Fq for the Hubbard 
model on the cubic lattice. Here the white region stands for 
values Fq > 10. No instability is found. 


instability was only found in the doped Mott insulator 
regime [l||. In that case, no intersite FM coupling is 
needed and the FM instability follows from Hund’s ex¬ 
change. 

On the contrary, an AF coupling suppresses the FM 
instability, and the value of Jo = 0.1 totally removes 
ferromagnetism as shown in Fig. [6] For the small values 
of U the Landau parameter Fq is negative in the entire 
regime of n, but then changes sign when U approaches 
U c . The exception here is the regime of low filling n < 
0.2, where Fq < 0 for U < U c . Larger AF exchange 
coupling Jq/U = 0.2 leaves only a narrow range of Fq < 0 
(not shown). 


Hubbard model. Here we also consider again the case of 
nearest neighbor interactions, 

_ f V for a bond (ij) with j £ M(i) , , 

y 0 otherwise. ’ ' ' 

where we use the same notation as in Eq. OOl). This case 
is particularly favorable for a charge instability near half 
filling, for instance in the form of a checkerboard state 
stable for large values of V as then the cost of any inter¬ 
site Coulomb energy oc V can be avoided. However, such 
a state does not occur as a consequence of the investi¬ 
gated instability of the uniform state. When V > 0, one 
finds even a stronger increase of Fq near half filling, and 
finally it becomes even larger than Fq = 10 in a broad 
range of filling n > 0.6 for the cubic lattice at V = 0.2 U, 
see Fig. [8] The uniform charge distribution is therefore 
more robust in the regime of n ~ 1, if V/U > 0. 

We have found that the uniform charge distribution is 
destabilized by attractive charge interactions V < 0, par¬ 
ticularly in the regime near quarter filling n ~ 0.5. At 
V = —0.217 the value of Fq decreases with increasing U 
for any filling n and this decrease is fastest near quarter 
filling. For U < U c one finds the charge instability at 
Fq = —1 in a broad range of n £ (0.045,0.93). This in¬ 
stability is related to the shape of the DOS and is easiest 
to realize at n ~ 0.42, there the DOS has a van Hove 
singularity. 

The data of Figs. [5] and ED suggest that in the case of 
charge response the regime near the metal-insulator tran¬ 
sition at half filling is robust and the Landau parameter 
Fq is here always enhanced, even for V <0. This find¬ 
ing confirms that electrons are strongly correlated and 
the system properties change radically at the metal-to- 
insulator transition at U c . Hence we inspect now the case 
n = 1 in more detail. We remark that the effect of finite 


B. Charge instability — parameter 

The symmetric Landau parameter Fq which stands for 
the charge response has to be evaluated numerically even 
at half filling, except for V = 0. As expected, Fq vanishes 
both for {7 = 0 and for n = 0, as Fq does. Otherwise, 
unlike the antisymmetric Landau parameter Fq which 
decreases with increasing U in the Hubbard model (at 
Vo = 0), the symmetric parameter Fq increases with U 
in the entire regime of filling 0 < n < 1. This increase is 
stronger near half filling, where Fq > 10 for U/U c > 0.7 
in a range of small doping away from half filling, see Fig. 
[3 At half filling the value of the positive Fq is given 
by Eq. (l39l) . It rapidly increases and finally diverges at 
the metal-to-insulator transition (we recall that for the 
simple cubic lattice U c ~ 16.04f). Away from n = 1 
the increase of Fq is moderate and it follows the same 
pattern as 1/z 2 in Fig. [2j being another manifestation of 
strong electron correlations near half filling. 

The increase of Fq with increasing U/U c is enhanced 
by a positive intersite Coulomb repulsion in the extended 



FIG. 8. (Color online) Landau parameter Fq for the extended 
Hubbard model on the cubic lattice. Increase of Fq with 
increasing U/U c is enhanced by positive V/U = 0.2; large 
values of Fq > 10 are found for n > 0.6. 
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moderate coupling: Fq is reduced for attractive V while 
it is enhanced for repulsive V, see Fig. [TO] The reduc¬ 
tion of Fq occurs only for sufficiently large — V and is 
visible in Fig. [TO] for V/U = —0.15, and beyond. As 
a result a minimum in Fq develops at U ~ 0.4f/ c , the 
minimal value of Fq decreases with increasing —V, and 
a charge instability may be found at the critical value 
V/U < —0.234, see the inset in Fig. [TO] The instability 
moves towards lower values of U with decreasing V when 
the minimum of Fq becomes deeper with further decreas¬ 
ing V. Particularly interesting is the non-monotonic be¬ 
havior of Fq with increasing U for V <0. We therefore 
suggest that a sufficiently strong intersite Coulomb at¬ 
traction —V/U > 0.234 is necessary to induce phase sep¬ 
aration. The instability is absent for repulsive V, where 
the uniform charge distribution is locally stable. 


FIG. 9. (Color online) Landau parameter Fq for the extended 
Hubbard model on the cubic lattice with attractive intersite 
interaction V/U = —0.2. Large values of Fq > 3 are found 
only near n = 1, while the charge instability Fq = — 1 oc¬ 
curs for a broad range of 0.045 < n < 0.93. Note that the 
instability line Fq = — 1 extends to n = 1~, and stops at 



FIG. 10. (Color online) Landau parameter Fq for the ex¬ 
tended Hubbard model on the cubic lattice at half filling 
(n = 1) for selected decreasing values of intersite Coulomb 
interaction V from top to bottom: V/U=0 (black line), 
V/U > 0 (blue) — V/U = 0.05 (solid line), V/U = 0.15 
(dotted) and V/U = 0.25 (dashed-dotted line), and V/U < 0 
(red) — V/U = —0.05 (solid line), V/U = —0.15 (dotted) and 
V/U = —0.25 (dashed-dotted line). The inset shows the in¬ 
stability value U inst jUc for V/U £ [—1-0, —0.2]. Its end point 
is marked by a solid circle. 


V first resembles somewhat that of finite J for weak to 


VI. DISCUSSION AND SUMMARY 

We have presented the consequences of finite intersite 
interactions: exchange J-ij and Coulomb V)j, and their 
impact on the instabilities of correlated electrons on a 
cubic lattice in the nondegenerate band described by the 
extended Hubbard model. While the Hubbard model was 
originally introduced, inter alia , to explain the metal¬ 
lic ferromagnetism, this idea fails as the model does not 
predict an instability of the metallic phase towards spin 
polarized (weakly ferromagnetic) states for the major¬ 
ity of lattices. Such an instability is expected in first 
place at half filling where electron correlation effects are 
strongest, so this case is of particular significance. How¬ 
ever, it is well known that the metal-to-insulator transi¬ 
tion happens typically for a lower value of U than the 
ferromagnetic instability, given by the divergence of spin 
susceptibility when the antisymmetric Landau parameter 
takes the value Fq = — 1. 

The analytic result derived in this paper shows that 
intersite exchange interactions change radically the re¬ 
sults for the metallic phase below the metal-to-insulator 
transition near U c . If these interactions are globally fer¬ 
romagnetic, a spin-polarized ground state occurs here be¬ 
fore the system undergoes the metal-to-insulator transi¬ 
tion, i. e., at U < U c . This result could be of importance 
to decide about the properties of some compounds with 
strongly correlated electrons, for instance, heavy fermion 
systems. 

Somewhat surprisingly we have found that the sym¬ 
metric (charge) Landau parameter Fq has a rather rich 
behavior when intersite Coulomb interactions are in¬ 
cluded. Small interactions do not change the value of Fq 
significantly. Most unexpected is the absence of a charge 
instability at half filling for moderate intersite Coulomb 
interactions, independently of the sign of the intersite 
Coulomb interaction, as it takes a sufficiently large at¬ 
tractive Coulomb interactions —V/U > 0.234 for a charge 
instability to occur. For smaller —V, it develops in the 
intermediate regime of doping, i. e., for 0.045 < n < 0.93. 
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It may be ex pec ted that the intersite Coulomb elements 
are repulsive JUlOjl , except for systems with strong pola- 
ronic effects, which gives no charge instability. Whether 
or not such interactions could be sufficiently attractive 
to cause an instability is an experimental challenge for 
future studies. 

Summarizing, we have shown that the Landau param¬ 
eters Fq and Fq are both sensitive to different intersite 
elements — Fq is modified only by exchange elements, 
while F 0 S only by Coulomb elements. The present the¬ 
ory predicts a ferromagnetic instability in a strongly cor¬ 
related metallic system with globally ferromagnetic ex¬ 
change. The derived analytic result for Fq Eq. © is 
very remarkable as it uncovers that for any lattice the 
Hubbard model at half Riling is at the verge of the fer¬ 
romagnetic instability, which is triggered by an infinites¬ 
imal ferromagnetic intersite exchange. This result pro¬ 
vides a new context for the original idea of Kanamori 
[I] who introduced the Hubbard model as the simplest 
model of itinerant ferromagnetism. 

At the same time, interesting behavior of the symmet¬ 
ric Landau parameter Fq was found for attractive near¬ 
est neighbor Coulomb interactions. In contrast to Fq, 
a more interesting result is found here away from half 
filling, where attractive interactions lead to charge insta¬ 
bilities signaling a tendency towards phase separation. 
We also remark that we investigated here only instabili¬ 
ties towards ferromagnetism or phase separation. Other 
instabilities may also occur at finite values of k which is 
an interesting subject for future research. 
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Appendix: Inverse propagator matrix 

The derivation of the inverse propagator matrix for 
lattices with inversion symmetry follows the one by Zim- 
mermann et al. [46|. To second order in the fluctuations 
the action S may be decomposed into two parts: 

5 = £(spin) +5 (charge) ; ( A .l) 

^charge) = £ ^„(- q), S£„(g)^„(g), (A.2) 

q,H,is 


with 1 < \i (is) < 6, and 

S {8pin) = ]T SM-Q)S S ^A<l)SM<l ), (A.3) 

-> v 


with 7 < ijl (y) < 12 including the three spin components. 
We use the following sequence of boson fields: 1/1 = e > 
2 = d', V>3 = d", 1 p 4 = Po, ^5 = Po, 1p6 = a, ip 7 = p x , 
V>8 = Px, = Py, IplO = Py, Ipll = Pz, 4>12 = Pz, after 
having introduced the short-hand notation d! = 5ft(d) and 
d” = 3(d). 

Regarding the spin sector at zero frequency we find: 


S 7 , 7 (k) = 5 9i9 ( k) = S 11A1 (k) 

1 

qT'O'-'fc) 


= a - P 0 + -Po Jg, 


£o£o 


d 2 z 

dpi 


hoX2 (k) 


[ dp z ’ ' 


Ss,s{k) = Sio,io(&) = S 12 ,i 2 (k) = -^Xo(k), 

S 7 ,a(k) = S 9 , 10 (k) = S n ,i 2 (k) 

= -Po-\zoXi(k)^- (A.4) 


Here we used the Fourier transform Eq. m of the 
intersite exchange elements. Notably, regarding th e sp in 
sector, the only difference to the Hubbard model |46| is 
the presence of Jj: in S 7 , 7 (k). 

Regarding the charge sector at zero frequency we ob¬ 
tain S c = S + S ^ with the non-vanishing matrix ele¬ 
ments of the symmetric matrix S: 

Si,i(k) = a + —Vo , 

*52,2 (k) = a — 2/3o + U — —Vo , 

S 3 , 3 (k) = a - 2 (3q + U — -Vo , 

*54,4 (k) = a — Po, 

Si, 6 (k) = e, 

S2,5(k) = —2d , 

S 2fi (k)=d, 

S4,5(k) = -Po , 

S^fiik) = po . (A.5) 


We introduced above S( z \k) = S(k ) + S^ v \k), with the 
non-vanishing matrix elements of S(k) defined as: 


Sy,v(k) = £qZo 


d 2 z 


S3,3(k) = £qzq 


dip y dip „ 
At, 1 / = 1,2,4 
d 2 z 


+ ( £ k ~ Q z lx2{k) 


dz dz 
dipy dip v 


dd” 2 


+ [ £ k-p z oX' 2 (k))—— 


dz dz* 


~ 1 Dz -* 

Sy,s{k) = --Z0 W Xi{k), 


where 








11 


At = 1,2,4 
S 5 ,5(k) = ~^Xo{k). 


(A.6) 


These elements are defined in terms of: 

jF{E %+p )-fF{E P ) 

Sp "H L k+pJ 

p,a 


Xn(k) = -^2{tf+t Up y 


Ek+p Ep 


x' 2 (k) = -J2(tp-h+p) 


2 M e Up)-M e p ) 

-%+p — Ep 


— y',tn-kf F ( E p)- 


"k Z_y “p 

P,(J 


(A.7) 


(£)=^ez 0 Vj:J^xi(>) 


®(fc)= 


cte . 


-d V-xoO) -2dz 0 —V^xi(k) + £y 


££? (*)=-£ dzoV^^-xi(k): 

sZ\k)=-l dV kXo(k), 


S$3(k)=-e v . 


(A.8) 


For fc = 0, the relation x' n (0) = (2i_F) n X'o(0) may be 
established |60| . 

Explicitly, using V _£ = Vg due to the inversion sym¬ 
metry of the lattice, the relevant matrix elements of the 
symmetric matrix read: 


sZ\k)=- 2 

s[ V 2 \k)=l 


c) z —* 

~e 2 V£xo(k) + 2ez 0 —V^xi{k) - £y 


edV~xo{k) +z 0 V m (k) ( e— -d— 


dz ,dz 


where we also introduced: 


£v = Vq n. 


(A.9) 


At zero frequency the field ifa decouples while for V = 0 
Zimmermann et aids result [461 ] is recovered. 
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